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Abstract 

A  propeller  lifting  surface  design  and  analysis  program  is  automatically  coupled  with  a  Euler 
/  Integrated  Boundary  Layer  Theory  (IBLT)  axisymmetric  flow  solver.  The  lifting  surface 
method  solves  the  localized  propeller  problem,  while  the  Euler /IBLT  solver  handles  the  global 
flowfield,  capturing  the  effective  inflow  problem.  For  viscous  flows,  the  boundary  layer  is 
constructed  based  upon  the  parameters  specified  by  the  IBLT  solution,  and  is  merged  with  the 
inviscid  Euler  flowfield.  The  robust  coupled  method  is  capable  of  representing  open  propellers, 
ducted  proulsors,  and  internal  flow  cases,  including  multi-blade  row  applications.  For  large 
axisymmetic  bodies,  the  user  may  specify  a  nominal  inflow,  and  the  coupled  method  is  used 
for  the  localized  propulsor  problem  only,  further  increasing  the  computational  efficiency.  The 
specified  nominal  flow  field  may  be  calculated  by  other  numerical  flow  solvers,  obtained  from 
experimental  results,  or  calculated  from  a  Euler/IBLT  solution  of  the  entire  body.  The  coupled 
code  is  an  extremly  rapid  flowfield  gridding,  calculation,  and  convergence  method,  which  allows 
an  order  of  magnitude  reduction  in  convergence  time  when  compared  to  the  current  efforts 
using  Reynolds  Averaged  Navier  Stokes  (RANS)  as  the  flow  solver.  Experimental  validation 
is  demonstrated  for  open,  ducted,  and  internal  flow  propulsors. 

Thesis  Supervisor:  Justin  E.  Kerwin 
Title:  Professor  of  Naval  Architecture 
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Chapter  1 


Introduction 

1.1  Overview 

The  goal  of  this  research  was  the  creation  of  a  fully  automatic  coupled  method  between  a 
lifting  surface  propeller  blade  design  and  analysis  program  and  an  Euler  throughflow  solver. 
A  coupled  technique  is  necessary  as  the  propeller  influences  the  flowfield,  and  the  flowfield 
influence  the  propeller.  Thus,  both  the  propeller  and  flow  field  calculation  must  be  iteratively 
converged  to  a  final  solution.  Currently,  coupled  design  and  analysis  of  marine  propulsors  is 
generally  accomplished  by  using  a  Reynolds  Average  Navier  -  Stokes  (RANS)  flow  solver,  the 
use  of  which  is  a  computational  intensive  task  requiring  a  high  degree  of  operator  knowledge 
and  experience.  A  Euler  solver  with  an  Integral  Boundary  Layer  Theory  (IBLT)  boundary 
layer  reconstruction  routine,  enables  solution  convergence  twenty  to  thirty  times  as  fast  as  the 
current  RANS  methods. 


1.2  Lifting-Surface  Propeller  Blade  Design  and  Analysis 

The  use  of  lifting  surface  theory  for  the  design  and  analysis  of  propellers  has  been  extensively 
documented  in  the  literature,  and  wall  be  only  summarized  here  for  the  sake  of  completeness[3],[ll],[12]. 
Propeller  Blade  Design  (PBD-14)  developed  by  Kerwin1  et  al  was  used  for  this  research.  A 
lattice  of  discrete  vortex  segments  is  placed  along  the  blade  mean  camber  surface,  the  hub,  the 

1  Justin  E.  Kerwin,  Profesor  of  Naval  Architecture,  Ocean  Engineering  Department,  MIT. 
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duct,  if  present,  and  the  trailing  wake  system.  Control  points  are  placed  a  particular  point 
within  each  grid  lattice.  Each  lattice  segment  has  an  appropriate  vorticity  strength,  such  that 
the  kinematic  boundary  conditions  are  satisfied  along  the  gridded  surface,  namely  the  total 
velocity  must  be  tangent  to  the  surface  at  each  control  point. 

In  propeller  design  and  analysis,  the  velocity  relation  is 

V  total  —  V  induced  V  effective  (1*1) 

where 


effective 


the  total  velocity  in  the  presence  of  an  operating  propeller 
the  velocity  induced  by  the  propeller  and  trailing  wake  vortex  distributions 
the  velocity  due  to  the  interaction  between  the  propeller  vortex  distributions 
and  the  presense  of  vorticity  in  the  inflow 


V  nominal  =  the  velocity  present  in  the  absence  of  the  propeller 


It  is  important  to  realize  that  neither  V  induced  nor  V  effective  are  physical  velocities  that 
may  be  measured.  Only  V  total  actually  exits  in  the  flow  field,  and  its  division  into  two 
components  is  done  solely  for  computational  reasons.  It  is  important  to  point  out  that  the 
above  listed  velocities  are  completely  independent  from  the  nominal  velocity  V  nominal  which 
is  the  velocity  present  in  the  absence  of  the  propeller.  Based  on  the  above  definitions,  it  is 
useful  to  note  that  if  no  vorticity  is  present  in  the  nominal  flowfield,  then  no  interaction  between 
the  propeller,  as  modeled  by  vortex  distributions,  and  the  inflow  are  necessary.  Thus,  in  this 
special  case ,  the  nominal  inflow  is  equal  to  the  effective  inflow.  This  is  never  the  case  in  reality 
since  viscous  effects  are  always  present,  however  it  often  serves  as  a  useful  check  of  propeller 
analysis  code  operation  when  run  in  an  inviscid  mode. 

The  lifting  surface  method  solves  the  propeller  problem  by  first  developing  an  influence 
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matrix  [JiVF]  which  gives  the  velocity  induced  by  each  vortex  lattice  segment,  with  an  assumed 
unit  strength,  upon  each  control  point.  When  this  influence  matrix  is  multiplied  by  the  actual 
vortex  segment  strength,  the  induced  velocity  at  every  control  point  is  known. 

[INF]  p]  =  Vinduced  (1.2) 

The  kinematic  boundary  condition  is  specified  by  dictating  that  the  total  velocity  normal 
to  the  grid  surface  must  be  zero.  Thus,  for  a  zero  thickness  case, 

[[inf]  [T]  T  V  e/Zectiuej  ^  —  0  (1-3) 

A  lifting  surface  code  may  be  used  for  either  a  design  of  a  new  propeller  or  an  analysis  of 
an  existing  propeller.  In  both  cases,  Equation  1.3  applies,  only  the  knowms  versus  unknowns 
vary.  In  the  case  of  propeller  design,  a  desired  radial  and  chordwise  loading  distribution  is 
prescribed,  and  the  blade  shape  is  manipulated  until  Equation  1.3  is  satisfied.  In  propeller 
analysis,  the  blade  shape  is  prescribed,  and  the  resulting  circulation  is  solved.  Blade  thickness 
effects  are  added  by  placing  discrete  source  lines  coincident  with  the  spanwise  blade  vortex 
lines.  The  resulting  propeller  forces  due  to  the  lifting  surface  and  thickness  are  calculated 
from  the  Kutta-Joukowski  and  Lagally  theorems,  respectively.  A  leading  edge  suction  force 
and  Lighthill  pressure  distribution  correction  are  applied  to  those  forces,  and  the  propeller’s 
sectional  viscous  drag  is  calculated  based  on  either  stripwise  two-dimensional  empirical  drag 
coefficients  or  a  stripwise  two-dimensional  integral  boundary  layer  calculation. 

1.3  Euler/IBLT  Throughflow  Solver 

MTFLOW  is  a  Euler/IBLT  throughflow  solver  developed  by  Drela2.  MTFLOW  incorporates 
a  streamline  curvature  method  is  a  powerful  method  for  solving  the  inviscid  flowfield  within  a 
channel  or  annular  passage.  As  the  method  is  inviscid,  a  boundary  layer  solution  must  be  used 
to  calculate  boundary  layer  quantities.  The  inviscid  streamlines  are  then  displaced  off  the  body 
by  the  displacement  thickness,  <!>*,  which  causes  a  redistribution  of  the  inviscid  streamlines. 

2 Mark  Drela.  Associate  Professor,  Department  of  Aeronautics  and  Astronautics,  MIT 
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Figure  1-1:  Finite  volume  geometry  for  throughflow  solver [2]. 

1.3.1  Euler  Throughflow  Solver 

MTFLOW  is  based  on  the  conservative  formulation  of  the  steady  state  Euler  equations  [1].  It 
relies  on  a  finite- volume  throughflow  method  based  on  the  streamline  curvature  method.  Figure 
1-1  displays  the  finit- volume  geometry  representation  used  within  MTFLOW. 

The  first  governing  equation  is  the  conservation  of  mass.  Thus,  across  adjacent  finite 
volumes  boundaries 

PxVmlSi  •  A  I  =  p2Vm2S2  *  A2  (1.4) 

The  conservation  of  momentum  is  written  as  the  steady  state  Euler  equation. 


P\  A  l  +  {Pi'OmiSi  •  Ai)vmiSi  +n  B  —  P2  -A  2  T  {p2Vm2S2  *  A2)vrn2s2  +  n+  B  +  (1.5) 

Conservation  of  energy  is  solved  for  in  terms  of  enthalpy 


■  -  1  px  +  2  ml 


ft  ,  1  2 
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Additionally  the  average  pressures  along  the  faces  pseudo-normal  to  and  along  the  stream- 


tube  must  be  equal. 


pa  +  p2  =  n+  +  ir 


MTFLOW  further  simplifies  equations  1.4  through  1.7.  The  meridional  flow  speed,  vm  is 
obtained  from  a  local  streamtube  conservation  of  mass. 


Um~  pA(2irr-BTe)  ^ 

The  streamwise  momentum  equation  solved  by  MTFLOW  is  of  the  same  form  as  Equation 


dP  +  pvmdvm  +  pvodvQ  +  Pd(AS)  —  pd(AW)  =  0  (1.9) 

The  change  in  work,  say  due  to  the  presence  of  a  propeller,  AW  enters  the  formulation 
through  the  Euler  turbine  equation. 


AW  =  J  Qd(rvo) 


(1.10) 


Q  =  blade  rotation  rate  ( - ) 

sec 


This  results  in  a  computationally  efficient  methodology  for  solving  the  inviscid  equations  of 
motion,  and  is  ideally  suited  to  couple  with  the  propeller  lifting  surface  code  PBD. 


1.3.2  Integral  Boundary  Layer  Theory  Equations 

Laminar  and  turbulent  boundary  layer  flows  are  governed  by  the  same  integral  momentum 
equation,  obtained  by  integrating  a  combination  of  the  continuity  and  simplified  cr-momentum 
equations  term  by  term  across  the  boundary  layer  [17]: 
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(1.11) 


where 


d£ 

dx 


Cf 

2 


fOO  ^  ^ 

6  —  momentum  thickness  —  /  —  (1  —  —  )dy 

J  o  Ue  Ue 

(1.12) 

*  f°°  u 

8  =  displacement  thickness  =  /  (1 - )dy 

JO  Ue 

(1.13) 

8* 

H  =  momentum  shape  factor  =  — 

V 

(1.14) 

du 

r w  =  wall  shear  stress  =  pz/7—  L 

dy 

(1.15) 

Cf  =  skin  friction  coefficient  =  ——r- 

pul 

(1.16) 

In  laminar  flow,  the  variables  6 ,  H:  and  Cf  can  be  reasonably  well  related  with  one- parameter 
velocity  profile  approximations.  For  marine  vehicles  moving  at  typical  speeds,  however,  the 
transition  from  laminar  to  turbulent  flow  within  the  boundary  occurs  near  the  bow,  and  the 
boundary  layer  remains  turbulent  over  the  remainder  of  the  vehicle.  The  turbulent-flow  pro¬ 
file  is  complicated  in  shape,  and  many  different  correlations  or  additional  relations  have  been 
proposed  to  effect  closure  of  Equation  1.11  [17].  The  closure  relations  used  within  MTFLOW 
are  those  derived  by  Swafford[15].  It  is  important  to  note  again  however,  that  MTFLOW  does 
not  solve  for  the  boundary  layer  velocity  profile,  but  rather  computes  the  boundary  layer  pa¬ 
rameters  and  displaces  the  inviscid  flow  field  off  the  body  by  S* ,  and  thus  computes  the  global 
flow  field  inviscidly. 

One  of  the  simplest,  and  thus  robust,  approaches  to  express  the  boundary  layer  velocity 
profile  and  recommended  for  general  use  by  White  [17],  is  the  ^ th  power  law: 


u_  _  /  y  y/7 

Ue  V<W 


(1.17) 
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where  £99,  the  boundary  layer  thickness,  is  the  value  of  y  at  which  the  boundary  layer 
velocity  is  equal  to  99  %  of  free  stream  velocity.  By  substituting  Equation  1.17  into  Equation 
1.13  and  replacing  the  upper  limit  of  infinity  of  the  integral  with  699,  the  relationship  between 
<5*  and  699  is: 


<$99  =  $8  (1.18) 

This  relation  is  useful  to  determine  the  boundary  layer  thickess  given  a  known  displacement 
thickness.  Then  Equation  1.17  may  be  used  to  determine  u  as  a  function  of  y . 

1.3.3  MTFLOW  Operation 

MTFLOW  consists  of  three  executable  programs;  MTSET,  MTFLO,  and  MTSOL.  MTSET 
creates  the  body  geometry  and  automatically  sets  up  the  initial  grid.  MTFLO  allows  the 
input  of  external  flow  field  information,  such  as  rotor  generated  swirl  or  additional  losses,  into 
the  flow  domain.  MTSOL  executes  the  flow  field  solution.  The  only  required  input  file  for 
MTFLOW  is  the  walls. xxx  file  which  is  used  by  MTSET  to  define  the  body  geometry.  The 
creation  of  the  walls. xxx  file  and  the  execution  of  the  MTFLOW  programs  are  sufficiently 
covered  in  Reference  [1]  and  will  thus  not  be  repeated  here.  Hereafter,  MTFLOW  refers  to  the 
entire  throughflow  solver,  while  MTSET,  MTFLO,  and  MTSOL  refer  to  the  individual  running 
of  the  respective  executable  programs. 


1.4  Objective 

This  thesis  provides  a  description  of  the  fully  automated  coupling  between  the  throughflow 
solver  MTFLOW  and  the  lifting  surface  propeller  blade  design  and  analysis  code  PBD.  An 
overview  of  the  process  is  first  presented,  followed  by  more  detailed  descriptions  of  the  coupling 
processs.  Finally,  validation  examples  are  presented. 


18 


Chapter  2 


Coupling  PBD  to  MTFLOW 

2.1  Process  Overview 

The  coupling  of  the  lifting  surface  propeller  blade  design  and  analysis  code  and  the  streamline 
curvature  code  follows  the  same  logic  as  utilized  by  Kerwin  et  al  [9]  when  coupling  the  propeller 
blade  design  and  analysis  code  with  a  RANS  code.  The  basis  of  Kerwin’s  coupling  was  the  use 
of  distributed  body  forces  within  the  RANS  domain  to  represent  the  presence  of  the  propeller. 
In  the  streamline  curvature  method,  the  propeller  is  modeled  as  a  swirl  introduction  device. 
Thus  added  swirl,  or  angular  momentum,  is  used  to  represent  the  presence  of  the  propeller 
in  the  flowfield.  This  technique  was  first  demonstrated  by  Renick  [14]  in  a  manual  coupling 
method,  and  the  present  research  is  an  extension  of  those  efforts  to  produce  a  fully  automated, 
robust  method  capable  of  solving  the  full  range  of  propulsor  types.  While  the  overall  concept 
and  code  execution  is  similar  to  that  used  by  Renick,  the  coupling  codes  have  been  completely 
rewritten  and  additional  changes  have  been  incorporated  into  the  flowfield  solver  to  affect  a 
correct  and  fully  automatic  coupling  technique. 

The  coupling  technique  is  started  by  solving  the  nominal  flowfield,  i.e.  the  body  with  no 
propeller  present.  For  viscous  flow,  the  boundary  layer  is  reconstructed  and  merged  to  the 
flowfield.  This  nominal  flow  field  in  then  passed  to  the  lifting  surface  code  which  solves  the 
localized  propeller  problem.  The  updated  circumferential  mean  induced  velocities  are  converted 
to  swirl,  and  are  then  inserted  into  the  global  flowfield  which  is  then  reconverged.  The  process 
is  then  repeated  until  both  the  flowfield,  and  lifting  surface  propeller  results  have  converged. 
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-Lifting  Surface 


BL2Body 

-Boundary  Layer 
Reconstruction 


- Inviscid flowfield  (to  8*) 
-Boundary  layer  parameters 


■Input  rVe  and  losses  into 
&  MTFLOW 


Figure  2-1:  Coupling  sequence  between  propeller  blade  design  code  (PBD)  and  the  throughflow 
solver  (MTFLOW). 


The  overall  sequence  is  shown  in  Figure  2-1. 


2.2  Running  the  Coupling 

The  overall  coupling  is  controlled  by  the  coupling  administration  file,  mtcouple  .inp.  The 
coupling  is  executed  by  two  computer  codes:  PBD2MT,  which  converts  the  PBD  output 
to  MTFLOW,  and  BL2BODY,  which  passes  the  updated  flowfield  to  PBD.  The  codes  are 
described  in  Table  2.1. 

The  coupling  administration  file,  mtcouple. inp,  contains  the  required  conversion  informa¬ 
tion  to  execute  the  coupling.  A  sample  administration  file  is  displayed  in  Table  2.2 

The  input  lines  have  the  following  uses  within  the  code: 

1.  REYNOLD’S  NUMBER:  Reynold’s  number  for  the  throughflow  domain.  An  inviscid 
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Code 

Description 

PBD2MT 

1)  Converts  the  PBD  output  circumferential  mean  induced  velocity  to  swirl. 

2)  Writes  out  tflow.xxx,  the  input  file  to  MTFLOW  containing  the 
induced  rV$. 

BL2BODY 

1)  Converts  the  MTFLOW  output  into  velocity  flowfield. 

2)  For  viscous  cases,  reconstructs  the  boundary  layer  profile. 

Table  2.1: 

Coupling  code  description 

LINE  1 

3 

1 

Reynolds  number 

LINE  2 

0.01 

! 

inlet  mach  number 

LINE  3 

1.0 

1 

Vship  used  in  PBD  to  calculate  Js 

LINE  4 

2.0 

1 

x  location  of  LE  tip 

LINE  5 

1.0 

r  location  of  LE  tip 

LINE  6 

XXX 

1 

MTFLOW  case  name 

LINE  7 

xxx . pbd 

| 

pbd  input  file  name 

LINE  8 

1 

I 

BLinput  toggle  (0=no , l=yes) 

LINE  9 

0 

| 

Nominal  velocity  toggle  (0=no,l=yes) 

LINE  10:  2 

1 

Number  of  blade  rows 

Table  2.2:  Example  coupling  administration  file 


case  is  forced  by  specifing  a  Reynold’s  number  less  than  10. 

2.  INLET  MACH  NUMBER:  For  compressible  flows,  the  Mach  number  is  defined  as  the 
inlet  velocity  divided  by  the  speed  of  sound.  For  water,  M  —  166Q.m.  .  For  general 

sec 

use  with  water  a  Mach  number  of  0.01  is  recommended.  This  prevents  round  off  errors 
internal  to  MTFLOW,  and  models  the  flow  as  virtually  incompressible. 

3.  VSHIP:  This  is  the  upstream  velocity  at  infinity.  For  coupling  purposes  it  is  the  fluid 
velocity  at  the  domain  inlet.  The  advance  coefficient,  J,  within  PBD  is  calculated  using 
Vship.  Traditionally  cases  are  non-dimensionalized  with  VshiP  =  1.0. 

4.  X  LOCATION  OF  LE  TIP:  The  axial  location  of  the  propeller  leading  edge  tip  in 
MTFLOW  coordinates. 

5.  R  LOCATION  OF  LE  TIP:  The  radial  location  of  the  propeller  leading  edge  tip  in 
MTFLOW  coordinates. 

6.  MTFLOW  CASENAME:  The  casename  of  the  walls  geometry  file  read  by  MTSET. 
The  generated  state  and  tflow.xxx  file  will  be  automatically  labeled  with  the  casename 
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extension. 


7.  PBD  INPUT  FILE  NAME:  Name  of  the  PBD  administration  file.  Used  to  read  the 
advance  coefficient  and  number  of  blades. 

8.  BLINPUT  TOGGLE:  Specifies  whether  an  inlet  boundary  layer  profile  is  present.  If  set 
to  1.  then  the  file  BLin.txt  is  read,  converted  to  entropy,  and  added  to  the  flowfield. 

9.  NOMINAL  VELOCITY  TOGGLE:  The  toggle  is  set  to  1  if  the  nominal  profile  is  required 
for  use  in  adjusting  the  specified  input  boundary  layer  to  match  a  desired  nominal  profile 
at  the  plane  of  the  propeller. 

10.  NUMBER  OF  BLADE  ROWS:  Specifies  the  total  number  of  rotor  and/or  stator  rows 
present  in  the  propulsor.  Thus,  a  value  of  1  describes  a  traditional  propeller,  and  a  value 
of  2  is  correct  for  a  pre-swirl/rot  or  combination.  Up  to  a  total  of  nine  blade  rows  are 
permitted. 
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Chapter  3 


PBD  to  MTFLOW  Conversions 

3.1  Program  Overview 

Information  is  passed  into  MTFLOW  by  means  of  the  tf  low .  xxx  file.  The  tflow  file  may  contain 
information  for  single  blade  rows,  multiple  blade  rows,  entropy  loss  due  to  the  presence  of  an 
upstream  boundary  layer,  or  any  combination  of  the  above.  This  link  of  the  coupling  therefore 
requires  the  transfer  of  the  necessary  information  into  the  tflow  file.  This  is  accomplished  by 
the  program  PBD2MT. 

3.2  Program  Operation 

3.2.1  Program  Flow  and  Required  Files 

After  the  successful  execution  of  both  MTFLOW  to  compute  the  nominal  flowfield  and  the 
initial  PBD  run,  PBD2MT  is  used  to  create  the  tflow. xxx  file.  As  successive  iterations  are 
performed,  grid  information  from  the  previous  MTFLOW  run,  and  the  updated  PBD  output, 
are  used  to  write  an  updated  tflow  file.  By  comparing  the  sizes  of  the  various  grid  files,  the 
program  automatically  determines  the  existence  of  an  open,  ducted,  or  internal  flow  propeller; 
and  incorporates  the  necessary  details  into  the  tflow  file.  The  required  input  files,  source, 
and  use,  are  explained  in  Table  3.1  The  program  operation  is  diagramed  in  Figure  3-1.  The 
details  of  more  complicated  situations  of  multiple  blade  rows  and  inlet  boundary  layer  profiles 
are  described  in  later  sections. 
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File 

Source 

Use 

xxx.pbd 

User 

Uses  propeller  advance  coefficient  to  compute  rotation  rate. 

mtcouple.inp 

User 

Contains  location  of  propeller  in  the  flow  domain  and  casename. 

PBDOUT.CMV 

PBD 

Contains  induced  cimcumfential  velocity. 

Gridxxxxx.dat 

MTSOL 

Contains  grid  information  from  previous  iteration  used  to 
compute  nondimensional  streamline  numbers. 

Table  3.1:  Input  files  required  by  the  program  PBD2MT. 


tflow.xxx  SWIRL.PBD 

SWIRL.MTF 


Figure  3-1:  Program  order  and  file  passing  when  running  PBD  to  MTFLOW.  The  end  product 
is  the  ascii  file  tflow.xxx,  which  is  the  required  input  for  MTFLO. 

3.2.2  The  tflow.xxx  file 

Information  is  passed  into  MTFLOW  by  use  of  the  tflow.xxx  file  which  is  read  by  MTFLO. 
Prescribed  distributions  of  tVq.  A S,  etc.,  are  specified  as  an  MxN  grid  of  points  arranged  in 
a  logical  rectangle  in  the  parameter  space.  In  the  current  coupling  routine,  only  vVq  is 
specified  in  the  file.  The  grid  is  input  as  a  set  of  “profile  blocks” .  The  grid  is  aligned  with  s 
the  component  perpendicular  to  the  flow,  and  t  the  component  parallel  with  the  flow.  Each 
block  has  a  constant  s  value,  and  thus  defines  quantities  as  a  profile  perpendicular  to  the  flow 
at  that  s  value  [1].  An  example  tflow  file  is  shown  in  Table  3.2. 
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Label 

#  t 

#  k 

z/L 

B 

r/L 

Omega  L/q 

rVt/qL 

Thk/L 

DS 

DH/a~2 

*  1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

Kstage(l) 

t(l,l) 

t(2,i) 

Nblade(l) 

X(l,l) 

X(2,l) 

Omega(l) 

Y(l,l) 

Y(2,l) 

RVT(1,1) 

RVT(2,1) 

THK (1,1) 
THKC2.1) 

DS(1,1) 

DS(2,1) 

DH(1,1) 

DHC2.1) 

t(N,l) 

X(N,1) 

Y(N,1) 

RVT(N,1) 

THKCN.l) 

DS(N,1) 

DH(N,1) 

Kstage(M) 

Nblade(M) 

Omega (M) 

t(N,M) 

X(N,M) 

Y(N,M) 

RVT (N , M) 

THK(N,M) 

DS (N , M) 

DH(N,M) 

Table  3.2:  Example  tflow.xxx  input  file  to  MTFLO 


The  contents  of  the  file  are  described  below.  Once  again,  items  11  and  12  are  not  determined 
in  the  current  coupling  technique,  but  are  features  which,  in  future  work,  may  be  utilized  for 
further  improvement.  The  reference  length,  Lref,  is  defined  as  1.0  in  the  MTFLOW  coordinate 
system.  The  points  are  specified  with  n  values  increasing  perpendicular  to  the  flow,  and  m 
values  increasing  parallel  to  the  flow. 

1.  LABEL:  Arbitrary  label  text  line. 

2.  #  ...:  Comment  label  lines. 

3.  *  1.0  Multiplier  line  which  is  applied  to  all  subsequent  values  in  the  associated  column. 
PBD2MT  sets  all  multipliers  at  1.0. 

4.  Kstage(m):  Number  of  the  rotor  or  stator  stage  which  is  described  by  the  profile  block. 
This  MTFLOW  function  is  not  used  in  the  coupling  method,  and  is  set  to  zero  by 
PBD2MT. 

5.  Nblade(m):  The  number  of  blades  present  in  the  current  s  line. 

6.  Omega(m):  The  rotation  rate  in  OLre//<fcn;. 
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7.  t(n,m):  The  normalized  streamfunction  £,  described  below. 

8.  X(n,m):  The  x  coordinate  defined  as  z/Lref. 

9.  Y(n,m):  The  y  coordinate  defined  as  r/Lref. 

10.  RVT(n,m):  The  swirl  defined  as  tVq / qiniLref  • 

11.  THK(n.m):  The  blade  thickness  Te/Lref. 

12.  DS(n.m):  The  added  entropy  AS  from  some  adiabatic  loss  process. 

13.  DH(n,m):  The  added  enthalpy  A H/a%nl  from  heat  release.  (  Not  applicable  to  propeller 
analysis.)  This  is  always  set  to  zero  by  PBD2MT. 

3.2.3  Calculation  of  the  Normalized  Streamfunction 

The  normalized  streamfunction,  t,  is  defined  as: 


t  =  tty) 


^  ~Wi 


(3.1) 


where,  in  the  MTFLOW  application,  ^  and  are  the  streamfunction  values  at  the  innermost 
and  outermost  streamsurfaces  of  the  flowpath  between  any  two  bodies,  or  elements,  in  the 


flowfield.  With  multi-element  configurations,  ipi and  t/>2  are  the  streamfunction  values  at  the 
bottom  and  top  of  one  flow  channel,  and  t  is  additionally  offset  by  the  channel  number,  so  that 


0  <  t  <1  between  the  first  and  second  element,  1  <  t  <  2  between  the  second  and  third  element, 
etc.  In  all  situations,  each  streamsurface  has  a  unique  t  value  associated  with  it  [1].  For  a 


given  z,r  location,  the  normalized  streamfunction  may  be  calculated  by: 


fo  p7rr2dr 

^  J  Jimin  r _ 

f max  cmr2dr 

^  ^min  ~ 

which  for  uniform  flow  reduces  to: 


(3.2) 


where 


(3.3) 
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Figure  3-2:  Streamline  tracing  required  to  find  nondimensional  streamline  number  for  any  given 
point  of  interest  in  the  flow  field. 


R  —  radius  of  point  of  interest 
Rmin  =  radius  of  bottom  of  local  flow  field 
Rma,x  —  radius  of  top  of  local  flow  field 

However,  uniform  flow  generally  only  exits  at  the  inlet  of  the  MTFLOW  domain.  Thus, 
the  following  simple  and  robust  approach  is  used  to  calculate  the  value  of  t  given  any  z,r  point 
of  interest.  The  flowfield  streamlines  from  the  previous  MTFLOW  iteration  are  available  in 
the  Gridxxx.dat  files.  The  applicable  Gridxxx.dat  file  is  selected,  and  the  streamline  above 
and  below  the  point  of  interest  is  determined.  The  fractional  difference  between  the  point  of 
interest  and  the  identified  streamlines,  is  traced  back  to  the  inlet,  where  uniform  flow  exists 
based  on  the  inlet  MTFLOW  boundary  conditions.  Figure  3-2  shows  the  utilized  streamline 
geometry.  Thus  the  following  equation  is  used  to  calculate  R : 
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B  =  0  B*0  B  =  0 


0=0  0*0  0=0 


Figure  3-3:  Field  parameter  s,t  grid  for  an  open  propeller  on  a  straight  shaft,  showing  doubled 
t  profiles  at  blade  leading  and  trailing  edge. 

This  calculated  R  value  is  then  used  in  Equation  3.3  to  determine  the  desired  streamfunction, 
t,  given  any  z,r  position. 

3.2.4  Tflow.xxx  Construction 

Figure  3-3  shows  a  general  s,£  grid  for  an  open  propeller  on  a  straight  shaft.  To  give  more 
flexibility  in  splining,  a  slope  break  or  discontinuity  in  either  tor  s  can  be  specified  by  two 
identical  t- value  or  s-  value  lines  in  each  block  of  the  tflow.xxx  file.  As  shown  in  the  example, 
the  blade  leading  and  trailing  edge  (5  =  2  and  5  =  4)  are  doubled  to  force  a  spline  break  in 
and  rV#.  Thus  two  identical  profile  blocks  with  identical  X,  Y  are  generated.  One  with 
and  rV$  set  to  zero,  and  the  other  with  and  rV$  set  to  the  actual  blade  leading 

(or  trailing)  edge  values.  Similarly,  a  slope  break  or  discontinuity  in  the  s  direction  may  be 
specified  by  two  identical  data  lines  within  a  given  profile  block.  This  is  done  at  the  blade  tip, 
to  enforce  a  zero  tVq  condition  outside  the  propeller. 
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Figure  3-4:  Sample  grid  demonstrating  the  relationship  between  the  various  parameters  within 
PBD2MT. 

3.3  Multi-Blade  Row  Method 

In  multiblade  row  applications,  the  induced  tVq  from  one  blade  row  is  convected  in  the  wake  to 
the  next  blade  row,  at  which  point  the  induced  tVq  from  the  next  blade  are  added  to  those  in 
the  wake  from  the  upstream  blade.  Kelvin’s  theorem  of  conservation  of  circulation  states  that 
for  any  ideal  fluid  the  circulation  is  constant  along  any  streamline.  Thus,  lines  of  constant  tVq 
are  convected  along  each  streamline.  Since  the  MTFLOW  grid  conforms  to  the  streamlines, 
this  coupling  method  is  ideally  suited  for  multiblade  row  applications.  Individual  tflow.xxx 
files  are  created  for  each  blade  row  in  their  own  subdirectory  as  described  in  Section  3.2.4.  Each 
tflow.xxx  is  then  moved  to  the  master  case  directory  and  renamed  based  on  the  order  of  the 
blade  rows.  For  exmple,  the  first  blade  row  tflow.xxx  file  is  renamed  tflow.Bl,  the  second 
blade  row  tflow.xxx  file  is  renamed  tflow.B2,  etc.  The  program  BUILDTFLOW  then  takes 
all  the  tflow.Bx  files,  determines  the  cumulative  tVq  at  each  s  line  location,  and  creates  the 
tflow.casename  file  containing  all  the  blade  row  information.  It  is  then  this  tflow.xxx  file 
which  is  read  by  the  program  MTFLO.  The  recommended  directory  structure  and  required  file 
organization  for  a  rotor/stator  case  is  displayed  in  Figure  3-5.  The  program  BUILDTFLOW 
can  handle  up  to  a  total  of  nine  blade  rows.  A  script  file  for  running  a  multi-blade  row  case  is 
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Figure  3-5:  Multiple  blade  row  coupling  directory  structure  and  program  flow  chart.  The 
process  is  best  controlled  by  a  batch  script  file. 

contained  in  Appendix  A. 

The  final  tflow.xxx  file  for  a  waterjet  is  shown  in  Figure  3-6.  Due  to  the  point  density  on 
the  blades,  the  multiple  points  appear  almost  solid.  Moving  left  to  right  through  the  flow  field, 
the  inlet  5-line  is  shown  first,  followed  by  the  closely  spaced  rotor  s-lines  containing  the  rotor 
generated  rVg.  The  mid-stage  5-lines  represent  the  wake  from  the  rotor.  The  second  set  of 
closely  spaced  5-lines  represent  the  stator,  with  their  associated  rVg  being  the  superposition  of 
the  rotor  wake  and  the  negative  swirl  induced  by  the  stator.  The  outlet  section  is  the  combined 
rotor/stator  wake,  and  as  the  stator  is  designed  to  remove  the  induced  swirl  of  the  rotor,  this 
final  wake  exhibits  only  a  small  amount  of  swirl. 
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Figure  3-6:  Screen  capture  of  MTFLO  displaying  the  tflow.xxx  file  information  for  a  ro¬ 
tor/stator  waterjet. 

3.4  Upstream  Boundary  Layer 

An  additional  feature  of  the  coupled  PBD/MTFLOW  method,  is  the  ability  to  model  only 
the  stern  portion  of  the  body  of  interest.  This  has  several  advantages,  some  of  which  are 
helpful  for  rapid  design  trade-off  studies,  and  some  of  which  are  required  to  successfully  model 
vehicles  with  a  significant  portion  of  the  propeller  within  a  boundary  layer.  The  idea  is  to 
only  grid,  and  spend  computational  effort,  on  the  stern  of  the  vehicle.  The  MTFLOW  domain 
inlet  is  located  forward  of  the  propeller  by  at  least  two  propeller  diameters.  Experimental  data 
indicates  that  the  presense  of  the  propeller  does  not  affect  the  flowfield  forward  of  this  point. 
Thus  this  size  flow  domain  is  the  minimum  required  to  adequately  capture  the  effective  wake [5]. 
Once  the  approriate  region  is  gridded  by  MTSET,  the  presence  of  a  boundary  layer  at  the  inlet 
to  the  domain  is  modeled  by  adding  an  entropy  loss  in  each  streamtube.  MTFLOW  is  then 
run  in  the  inviscid  mode,  but  with  the  boundary  layer  captured  as  an  entropy  loss.  While 
the  boundary  layer  is  assumed  to  convect  and  evolve  appropriately  as  the  shape  of  the  body 
changes,  the  additional  viscous  drag  effects  (which  would  normally  thicken  the  boundary  layer) 
occurring  between  the  inlet  plane  and  the  propeller  plane  are  not  included  in  the  inviscid  mode. 
However,  this  effect  is  considered  minimal  in  comparison  to  the  overall  accuracy  and  speed  of 
this  coupled  method,  especially  given  the  fact  that  the  action  of  the  rotating  blade  rows  will 
tend  to  accelerate  the  fluid  and  cancel  much  of  the  additional  boundary  layer  thickness  due  to 
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friction. 

Once  the  nominal  flow  is  calculated,  the  propulsor  induced  swirl  may  be  added,  and  the 
total  flow  field/propeller  problem  solved  iteratively.  The  speed  of  this  “stern  only”  method 
permits  a  rapid  analysis  of  multiple  propulsor  configurations  for  design  space  exploration.  The 
other  significant  advantage  of  this  feature,  is  that  MTFLOW  can  now  model  the  contribution  to 
the  flow  field  by  the  portion  of  the  propeller  within  the  boundary  layer  displacement  thickness. 
If  run  as  designed,  a  significant  portion  of  the  propulsor  on  a  submerged  body  would  be  within 
the  displacement  thickness,  which,  by  definition  is  outside  the  MTFLOW  domain.  The  desired 
nominal  profile  may  be  provided  by  other  flow  solvers,  such  as  RANS,  experimental  data,  if 
available,  or  MTFLOW  itself  may  be  used  to  model  the  full  body  and  run  in  the  viscous  mode. 


3.4.1  Velocity  to  Entropy  Conversion 

To  convert  a  velocity  profile  into  an  entropy  distribution,  we  start  with  the  entropy  equation: 


S  =  In 


Po 


For  small  changes  (low  Mach  limit), 


(3.5) 


Using  the  relation, 


2P  Po 


(3.6) 


pU2  =  7 pM2 


(3.7) 


and,  for  the  case  of  low  Mach  numbers  ( i .  e,  water)  assume  that  the  pressure  is  approximately 
equal  to  the  absolute  total  pressure 


so  that 


P~Po 


(3.8) 
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Inlet  Boundary  Layer  Profile  Specified  as 


0.5  0.6  0.7  0.8  0.9  1 


Vx 

Figure  3-7:  Straight  shaft  with  jth  power  law  boundary  layer  specified  by  an  entropy  loss. 


(3.9) 


While  the  use  of  the  specific  heat  ratio,  7,  and  Mach  number,  M,  with  water  flow  are 
admittedly  foreign,  their  use  is  important  to  correctly  couple  with  MTFLOW.  A  specific  heat 
ratio  corresponding  to  that  of  air,  7  =  1.4,  is  required  to  ensure  consistency  with  the  built  in 
MTFLOW  methodology.  The  Mach  number  is  as  specified  in  the  coupling  administration  file, 
mtcouple . inp. 

To  validate  the  use  of  entropy  as  a  boundary  layer  modeling  tool,  a  straight  shaft  test  case 
was  conducted  with  a  jth  power  law  boundary  layer  present.  The  results  are  demonstrated  in 
Figure  3-7.  The  agreement  between  the  desired  boundary  layer,  and  that  convected  through 
MTFLOW,  displays  the  effectiveness  of  simulating  a  boundary  layer  as  an  entropy  loss. 
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3.4.2  Boundary  Layer  Entropy  Modeling  on  Submerged  Body 

When  the  same  procedure  is  conducted  on  a  body  of  revolution,  the  boundary  layer  profile  at 
the  domain  inlet  matches  the  specified  profile,  but  the  boundary  layer  downstream  develops  dif¬ 
ferently  than  the  experimentally  measured  profile.  This  is  presumably  a  result  of  the  increasing 
invalidity  of.the  pressure  assumption  of  Equation  3.8  as  the  curvature  of  the  body  increases. 
However,  to  ensure  that  only  velocity  data  is  required  for  the  case  setup,  it  is  not  desirable 
to  introduce  a  pressure  correction.  To  overcome  this  deficiency,  a  procedure  is  available  to 
tune  the  specified  inlet  velocity  profile,  and  the  associated  added  entropy  loss,  to  result  in  the 
correct  nominal  flow  field  in  the  propeller  plane.  Thus,  while  the  resultant  inlet  flow  domain 
may  differ  from  the  measured  profile  there,  as  the  flow  develops  the  correct  nominal  profile  is 
present  in  the  region  of  the  propeller.  Once  the  necessary  inlet  profile  is  determined,  it  remains 
constant  for  use  with  the  propeller  present. 

The  procedure  and  program  used  to  scale  the  inlet  boundary  layer  to  achieve  the  desired 
nominal  profile  at  the  propeller  is  displayed  in  Figure  3-8.  Haung  Body  1  (DTMB  Model 
5225-1)  was  used  to  validate  this  approach  [8].  The  file  format  of  the  boundary  layer  profiles, 
BLin .  txt  at  the  domain  inlet  and  Exp977nom .  dat  at  the  propeller  plane,  are  shown  in  Table 
3.3.  MTFLOW  wall  correctly  determine  the  radial  velocity  component  based  on  the  body 
shape,  so  only  the  axial  velocity  component  must  be  specified  in  the  boundary  layer  profile 
files.  The  boundary  layer  A S  loss,  is  added  in  the  same  manner  as  the  blade  induced  tVq  and 
appears  as  an  MTFLO  input  variable  DS  as  shown  in  Table  3.2.  The  BLINPUT  TOGGLE  in 
the  coupling  administration  file  mt couple .  inp  must  be  set  to  1  to  include  the  A S  loss.  During 
the  setting  of  the  inlet  profile  to  match  the  desired  nominal  profile  at  the  propeller  plane,  a 
PBDOUT.CMV  file  is  required  to  specify  the  blade  location;  however,  no  rVe  is  desired.  This  too 
is  controlled  by  the  mtcouple.inp  file  by  setting  the  NOMINAL  VELOCITY  TOGGLE  to  1 
during  this  scaling.  To  prevent  overshooting  the  desired  profile,  a  damping  coefficient  of  0.25 
is  included  in  each  scaling  step.  The  scaling  is  accomplished  for  each  streamline  as  follows: 


VBLini'lp)  new  —  Vb Unsold  X 


1  ^VmTFLOw{ i>)  ~  YoesiredW 
4  \  VmtflowO’) 


where 


(3.10) 
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NumPts 

R(l) 

V(l) 

R(2) 

V(2) 

R(3) 

V(3) 

R(NumPts) 

V(NumPts) 

Table  3.3:  Required  format  of  boundary  layer  file.  The  number  of  points  in  the  boundary  layer 
are  listed  in  the  first  line,  followed  by  each  R,  V  point,  where  V  is  the  axial  velocity  component. 


VsLin{ VO  new 
VBLinWold 
VmtflowM 
^Desiredi^P) 


=  updated  streamline  axial  velocity  in  boundary  layer  input  file  BLin.txt 

—  previous  streamline  axial  velocity  in  boundary  layer  input  file  BLin.txt 

—  streamline  axial  velocity  at  propeller  plane  in  MTFLOW 

=  streamline  desired  axial  velocity  at  propeller  plane  in  file  Exp977.nom 


Validation  of  this  method  was  conducted  with  Huang  Body  1  (DTMB  Model  5225-1).  Fig¬ 
ure  3-10  shows  the  required  inlet  profile  which  results  in  the  current  boundary  layer  in  the 
propeller  plane  shown  in  Figure  3-9.  The  complete  stern  section  of  the  body  and  boundary 
layer  growth  is  shown  in  Figure  3-11.  Of  note,  a  far  field  setting  of  a  constant  pressure  jet 
boundary  in  MTSOL  was  required  for  stable  convergence.  This  is  imposed  by  setting  the  far 
field  type  to  “3”  in  the  solution  parameter  settings. 

3.4.3  Gridding  Requirements  for  Stern  Section  Modeling 

The  automatic  gridding  procedure  contained  in  MTSET  was  designed  for  wings,  ducts,  or 
internal  flow.  Consequently,  when  directly  applied  to  the  stern  section  of  a  submerged  body, 
the  grid  smoothing  routines  may  have  difficultly  in  the  inlet  region  adjacent  to  the  body.  The 
bottom  left  corner  may  drift  above  the  next  streamline  as  shown  in  the  top  of  Figure  3-12.  By 
introducing  a  spline  break  as  a  doubly  specified  body  point  as  demonstrated  in  Reference  [1], 
and  artificially  lowering  the  inlet  body  point  slightly,  correct  gridding  is  accomplished.  This 
ensures  the  bottom  streamtube  has  positive  area,  and  thus  an  initial  positive  velocity.  The 
initial  grid  is  viewed  in  the  TECPLOT  formatted  file  ORGGRID  .tec.  The  existence  of  a  negative 
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mtcouple.inp 
PBDOUT.CMV  * 


**  Created  using  TecPlot polyline  extraction 


Figure  3-8:  Program  operation  and  file  passing  to  scale  the  specified  entropy  loss  to  match  the 
desired  nominal  profile  at  the  propeller  inlet  plane. 
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Figure  3-10:  Resulting  nominal  velocity  profile  at  X/L  =  0.914  in  MTFLOW  with  the  required 
entropy  loss  to  correctly  match  the  nominal  profile  at  the  propeller  location  of  X/L  =  0.977 
compared  to  experimental  profile  at  X/L  =  0.914. 


Inlet  Body  Points  Entered  at  Actual  Location 


Inlet  Body  Points  Adjusted  for  Correct  G ridding 


20.06  20.07  20.08  20.09  20.1 


Figure  3-12:  Grid  comparsion  showing  drift  of  bottom  left  corner  (top)  when  inlet  body  points 
entered  at  actual  location  versus  adjusted  body  points  (bottom). 

area  can  at  a  minimum  result  in  localized  flow  difficultly,  and  in  the  extreme,  can  lead  to  a 
negative  temperature  and  the  subsequent  crashing  of  the  program  MTSOL.  Fortunately,  once 
a  successful  grid  is  created,  it  is  sufficient  for  all  subsequent  coupled  runs.  A  corrected  grid  is 
shown  in  the  bottom  half  of  Figure  3-12.  The  localized  point  adjustment  has  no  impact  to  the 
overall  flow  downstream  in  the  propeller  region. 
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Chapter  4 


MTFLOW  to  PBD  Conversions 


4.1  Program  Overview 

The  link  from  MTFLOW  to  PBD  is  accomplished  by  the  program  BL2BODY.  The  output 
from  BL2B0DY  is  the  Tecplot®  formatted  flow  field  velocity  file,  VEL JOIN. tec,  which  is 
read  by  VELCON  [4],  the  traditional  velocity  conversion  program  used  to  create  the  updated 
flow  field  for  use  by  PBD.  For  an  inviscid  case,  BL2BODY  simply  converts  the  MTFLOW 
output  to  the  correct  format  required  by  VELCON.  In  a  viscid  case,  the  MTFLOW  output 
only  contains  inviscid  velocities  for  a  body  offset  by  the  displacement  thickness,  and  thus  the 
flowfield  must  be  expanded  to  the  body  boundary,  and  an  appropriate  boundary  layer  velocity 
profile  reconstructed. 


4.2  Program  Operation 

The  overall  program  execution  order  is  displayed  in  Figure  4-1.  The  coupling  administration 
file,  mtcouple .  inp,  is  read  to  determine  whether  a  viscid  or  inviscid  case  exits.  If  the  specified 
Reynolds  number  is  less  than  ten,  an  inviscid  case  is  assumed.  For  the  inviscid  case,  the 
VEL  JOIN,  tec  file  is  written  with  the  velocities  referenced  to  free  stream  velocity,  based  on  the 
velocity  and  Mach  number  as  referenced  in  the  coupling  administration  file.  For  a  viscid  case, 
a  boundary  layer  reconstruction  occurs.  The  beginning  and  ending  points  of  the  bodies,  and 
the  number  of  total  bodies  in  the  flow  field,  are  determined  from  the  walls  .xxx  file.  The  body 
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VELJOIN.tec 

l 


VELCON 

V _ ✓ 

l 

restart.vel 

Figure  4-1:  Program  order  and  file  passing  when  running  MTFLOW  to  PBD.  The  end  product 
is  the  TECPLOT  format  file  restart.vel  (or  other  name  as  specified  during  VELCON  execution), 
which  is  the  required  input  for  PBD. 

geometry  is  read  from  the  ORRIGGRID.tec  file.  As  this  file  is  created  by  MTSET,  before  a  flow 
solution  is  created,  this  grid  conforms  to  the  actual  body  location.  The  original  work  by  Renick 
[14]  utilized  a  Swafford  boundary  layer  profile,  however,  the  present  coupling  reverted  to  a  pure 
jth  power  law  reconstruction.  While  this  is  less  accurate  in  high  pressure  gradient  regions, 
it  is  more  robust,  and  was  thus  more  consistent  with  the  goal  of  a  rapid,  stable,  automatic 
coupling  method.  Certainly,  the  addition  of  a  more  advanced,  yet  robust  boundary  layer 
reconstruction  technique  would  improve  the  coupling  accuracy.  Of  note,  the  recommended,  and 
often  necessary  method  of  characterizing  a  boundary  layer  due  to  the  presence  of  a  significant 
displacement  thickness  in  the  propeller  region  for  a  viscid  case,  is  the  entropy  loss  method 
described  in  Section  3.4.2.  When  a  boundary  layer  is  created  using  an  entropy  loss,  MTFLOW 
is  run  in  the  inviscid  mode,  and  the  drawbacks  of  the  ^ th  power  law  are  irrelevant. 

The  velocity  output  from  MTSOL  is  written  in  the  file  OUTVEL.tec,  and  for  a  viscous  case, 
the  boundary  layer  information  is  written  in  the  file  OUTBL .  tec.  A  sample  OUTVEL .  tec  velocity 
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field  is  shown  in  Figure  4-2.  The  maximum  displacement  thickness,  max  5*,  is  determined  from 
the  OUTBL.tec  file,  and  following  the  derivation  in  Section  1.3.2,  the  maximum  boundary  layer 
thickness,  max  <$99,  is  calculated  by  Equation  1.18.  The  streamline  immediately  outside  the 
max  699  is  considered  the  lowest  fully  inviscid  streamline.  The  MTFLOW  grid  is  cut  at  this 
location,  and  boundary  layer  reconstruction  is  conducted  between  the  body  and  the  lowest  fully 
inviscid  streamline.  First,  new  high  density  grid  lines  are  faired  in  between  the  lowest  fully 
inviscid  streamline  and  the  body.  Then,  velocities  are  assigned  to  the  new  grid  nodes.  If  an 
individual  node  is  above  the  local  699,  then  the  OUTVEL.tec  grid  is  interrogated  to  determine 
the  local  velocities.  If  an  individual  node  is  between  the  body  and  the  local  £99  then  the  local 
axial  and  radial  velocity  components  are  calculated  in  accordance  with  the  ^th  power  law  of 
Equation  1.17.  Between  the  local  <5*  and  local  699  the  tangential  velocity  is  assigned  the  value 
in  the  MTFLOW  domain  as  contained  in  OUTVEL.tec.  Between  the  body  and  the  local  <5* 
no  tangential  velocity  data  exists  within  MTFLOW.  To  resonable  approximate  the  tangential 
velocity  in  this  region,  the  tangential  velocity  is  linearly  interpolated  between  the  actual  value 
at  the  edge  of  the  MTFLOW  domain,  <5*,  and  an  assumed  value  at  the  wall  of  50  %  of  that  at 
the  edge.  Figure  4-3  shows  an  example  VELJOIN.tec  file  written  by  BL2BODY  with  the  new 
viscid,  dense  grid  joined  to  the  MTSOL  calculated  fully  inviscid  grid  region. 

To  compare  the  boundary  layer  reconstruction  with  the  experimental  data,  velocity  profiles 
were  extracted  from  the  converged  viscous  Huang  Body  1  nominal  profile  case.  Figure  4-4 
shows  the  results  at  four  X/L  locations.  As  expected,  the  ~th  power  law  fails  to  accurately 
represent  the  experimental  profile  in  the  tapered  stern  section.  However,  it  still  serves  as  a 
useful  tool  in  the  absence  of  more  detailed  information,  and  for  use  along  straight  shafts,  ducts, 
and  internal  flow  cases. 
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freestream  velocity 


Figure  4-2:  Sample  viscous  flowfield  output  from  MTFLOW. 


Figure  4-3:  Sample  BL2BODY  grid  with  boundary  layer  reconstruction  grid  added  to  the 
inviscid  MTFLOW  grid. 
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1.1 


X/L  =  0.755 


1.1 


X/L  =  0.846 


Figure  4-4:  Comparison  of  reconstructed  boundary  layer  profile  and  experimental  measurements 
for  Huang  Body  1.  As  expected,  the  accuracy  of  the  ~ th  power  law  diminishes  in  the  highly 
tapered  stern  region. 


Chapter  5 


Validation 


Code  validation  was  conducted  on  the  following  cases: 

1.  Open  Propeller 

2.  Ducted  Propeller 

3.  Water  jet 

4.  Submerged  Body 


The  open  water  comparison  used  propeller  4119  in  comparison  with  the  1998  International 
Towing  Tank  Committee  (ITTC)  tests[13] .  The  ducted  propeller  KA-455  from  the  Netherland 
Ship  Model  Basin  (NSMB),  Kaplan  series  ducted  propulsor  [16],  was  utilized  for  the  ducted 
case.  Internal  flow  was  validated  by  the  Water  Jet-21,  as  tested  in  the  Marine  Hydrodynamics 
Laboratory  (MHL)  at  the  Massachusetts  Institute  of  Technology  (MIT)  [10].  Huang  Body  1 
(DTMB  Model  5225-1)  with  propeller  4577  [5]  was  used  for  the  submerged  axisymmetric  body 
comparison. 


5.1  Open  Propeller 

Propeller  4119  on  a  straight  shaft  was  used  as  the  initial  validation  of  the  coupling  method. 
Additionally,  the  final  results  were  compared  with  the  experimental  results  from  the  1998  ITTC 
tests  as  the  open  propeller  test  case  [13].  The  inviscid  axial  effective  velocity  was  calculated 
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and  is  show  in  Figure  5-1.  Given  a  nominal  axial  inflow  of  1.0,  the  axial  effective  inflow 
velocity  should  be  equal  to  1.0  on  the  entire  blade.  The  near  unity  results  demonstrate  the 
effectiveness  of  the  coupling  method.  An  additional  useful  coupling  check  is  to  compare  the 
tangential  trailing  edge  velocities  between  PBD  and  MTFLOW.  This  is  synonymous  with 
comparing  the  circulation  distributions.  Thus,  the  trailing  edge  circumferential  mean  induced 
tangential  velocities  output  from  PBD  are  compared  to  the  tangential  velocity  at  the  same 
spacial  location  in  the  MTFLOW  calculated  flow  domain.  This  is  shown  in  Figure  5-2.  The 
exact  agreement  between  both  domains  is  an  important  verification  that  the  coupling  is  being 
carried  out  correctly.  Propeller  4119  results  are  compared  with  the  1998  ITTC  experiment 
results  in  Figure  5-3.  MTFLOW  was  run  in  the  viscous  mode,  allowing  boundary  layer  growth 
on  the  shaft.  PBD  was  run  with  both  the  sectional  drag  coefficients  and  thickness  effects 
included.  A  relatively  dense  lifting-surface  grid  of  35  x  35  vortices  was  required  for  grid 
convergence1 .  As  expected,  the  accuracy  of  the  solution  diverges  at  off  design  J  values  due  to 
the  inability  of  the  lifting  surface  method  to  adequately  handle  extreme  blade  angles-of- attack. 


5.2  Ducted  Propeller 

A  ducted  propulsor  validation  was  conducted  using  the  Kaplan  KA-455  with  nozzle  20  as  tested 
in  the  Netherlands  Ship  Model  Basin[16].  For  comparative  purposes,  the  thrust  value  contains 
only  the  rotor  generated  thrust,  and  does  not  contain  the  nozzle  effect.  Figure  5-4  is  the 
KA-455  as  modeled  in  the  lifting  surface  tool.  MTFLOW  was  run  for  both  an  inviscid  and 
viscid  case.  Once  again,  the  PBD  analysis  included  both  the  sectional  drag  coefficients  and 
thickness  effects.  The  results  are  shown  in  Table  5.1.  The  results  demonstrate  the  accuracy 
of  the  method  and  the  minimal  effect  that  the  viscous  boundary  layer  on  the  shaft  and  duct 
have  on  the  overall  solution. 

^rid  convergence  means  that  a  further  increase  in  the  grid  density  does  not  affect  the  results. 
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P4119  Inviscid  Axial  Effective  Velocity 
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Figure  5-1:  Countor  plot  of  the  effective  axial  velocity  in  the  presense  of  a  uniform  nominal 
inflow  velocity  of  1.0  for  P4119  on  an  inviscid  straight  shaft. 


4119  Trailing  Edge  Vtheta  Comparison  between 


Figure  5-2:  Comparison  of  the  circumferential  mean  tangential  velocity  at  the  blade  trailing 
edge  of  P4119  from  MTFLOW  and  PBD. 
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4119  Coupled  MTFLOW  vs  1998  ITTC  Experiment 


Figure  5-3:  Comparison  of  P4119  MTFLOW/PBD  results  with  1998  ITTC  Experiment.  PBD 
was  run  with  both  sectional  drag  coefficients  and  thickness  effects  included  and  with  a  vorticy 
grid  lattice  of  35  x  35. 


Figure  5-4:  Kaplan  KA-455  ducted  propulsor  as  modeled  in  PBD. 


MTFLOW/PBD 

MTFLOW/PBD 

Experimental 

Re 

inviscid 

500,000 

kt 

0.245 

10K  Q 

0.397 

0.392 

0.400 

Table  5.1:  Comparison  of  KA-455  inviscid  and  viscid  MTFLOW/PBD  rotor  results  with  ex¬ 
perimental  data  at  design  J=0.36. 
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Figure  5-5:  Contour  plot  of  circumfential  mean  tangential  velocity  of  WaterJet  WJ21  with 
streamlines  superimposed.  Flow  is  from  left  to  right.  The  rotor  leading  edge  is  located  at 
approximatly  X  =  0.2  and  the  stator  leading  edge  is  located  at  approximatly  X  =  1.0. 

5.3  WaterJet 

Internal  flow  verification  utilized  the  Waterjet  WJ21  tested  in  the  Marine  Hydrodynamics 
Laboratory  at  the  Massachusetts  Institute  of  Technology [10].  The  torque  and  thrust  of  the 
rotor  were  compared  to  an  equivalent  RANS  calculation  and  with  the  experimentally  measured 
torque.  As  with  the  ducted  case,  the  thrust  of  the  rotor  does  not  correspond  to  the  net  thrust 
produced  by  the  propulsion  system,  but  does  serve  as  a  useful  quantitative  comparison.  A  plot 
of  the  circumferential  mean  tangential  velocities  with  streamlines  superimposed  is  displayed  in 
Figure  5-5.  The  results,  including  a  computational  time  comparison,  are  shown  in  Table  5.2. 
This  particular  case  demonstrates  the  relative  usefulness  of  the  MTFLOW/PBD  method.  The 
results  are  quite  close,  in  spite  of  the  fact  that  this  is  a  particularly  complex  flow  case.  The 
reduction  in  computing  time  from  12  hours  to  30  minutes  is  very  substantial.  As  the  flow  in 
the  test  section  of  the  waterjet  is  dominated  by  potential  flow,  this  case  was  run  in  the  inviscid 
mode.  For  design  applications,  the  boundary  layer  growth  of  the  upstream  flowpath  would  be 
perfectly  suited  to  be  modeled  via  an  entropy  loss.  This  would  permit  the  blade  design  to  be 
conducted  within  the  actual  flow  field  expected,  while  not  requiring  the  complete  gridding  of 
the  piping  system  as  is  currently  required  by  RANS  methods. 
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MTFLOW/PBD 

RANS/PBD 

Experimental 

kt 

1.92 

1.87 

K  Q 

0.409 

0.396 

0.424 

Time 

30  min 

12  hours 

Table  5.2:  Comparison  of  waterjet  WJ21  MTFLOW/PBD  and  RANS/PBD  results  with  ex¬ 
perimental  data. 


5  10  x/R  15  20 

Figure  5-6:  Representation  of  Huang  Body  1  (DTMB  Model  5225-1). 

5.4  Submerged  Body 

This  test  case  is  one  of  a  series  of  axisymmetric  bodies  tested  by  the  David  Taylor  Model  Basin 
by  Huang  et  al  [7],  [8],  [6],  [5].  The  same  forebody  (DTMB  Model  5225)  was  used  for  all 
experiments  while  a  series  of  afterbodies  with  an  increasing  degree  of  taper.  These  bodies  have 
been  used  extensively  for  validating  solutions  to  the  effective  wake  problem  using  analytic  and 
numerical  methods.  The  afterbody  considered  here,  Afterbody  1,  is  a  non-separating  stern 
with  a  low  tailcone  angle.  A  profile  view  of  the  experimental  body  (DTMB  Model  5225-1)  is 
shown  is  Figure  5-6.  It  was  tested  in  the  presence  of  an  open  rotor  in  wind  tunnel  and  towing 
tank  facilities.  An  existing  seven-bladed  propeller  (P4577),  with  a  diameter  of  54.5  %  of  the 
hull  diameter  is  mounted  at  x/L  =  0.983.  As  the  actual  body  boundary  layer  displacement 
thickness  includes  a  significant  portion  of  the  propeller,  the  entropy  loss  method  of  Section  3.4.2 
is  used  in  this  application. 

To  correctly  set  the  nominal  profile,  experimental  boundary  layer  profiles  were  selected  at 
two  locations.  The  first,  located  at  X/L  =  0.914,  served  as  the  initial  inlet  plane  profile,  and 
the  second,  located  at  the  propeller  inlet  of  X/L  —  0.977,  served  at  the  goal  nominal  profile. 
The  inlet  plane  profile  was  then  adjusted  in  accordance  with  the  procedure  of  Section  3.4.2  until 
the  MTFLOW  nominal  propeller  plane  inlet  profile  matched  the  experimental  profile.  Once  set, 
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Nominal  and  Total  Velocity  Comparison  of 

MTFLOW  and  Experiment 

Huang  Body  1  with  P4577,  J=  1 .25,  at  X/L  =  0.977 
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Figure  5-7:  Huang  Body  1  nominal  and  effective  velocity  comparison  of  MTFLOW  and  exper¬ 
iment  at  X/L  =  0.977  with  inlet  entropy  loss  specified  at  X/L  =  0.914.  The  propeller  leading 
edge  tip  is  located  at  R  =  0.545. 

the  inlet  profile  remained  fixed  during  the  propeller  analysis.  The  nominal  and  total  propeller 
plane  profiles  are  shown  in  Figure  5-7.  The  calculated  and  experimental  thrust  and  torque 
coefficients  are  contained  in  Table  5.3.  Of  note,  the  listed  experimental  torque  coefficient  is 
not  purely  a  measured  value.  The  propulsor  used  in  the  self-propulsion  experiment  was  a 
stock  propulsor  and  was  not  operated  at  its  design  angle  of  attack.  The  experimental  Kq  was 
calculated  by  a  propulsor  performance  prediction  computer  program  based  on  hydrodynamic 
pitch  angles  tan  which  were  iteratively  scaled  based  on  the  ratio  of  measured  to  predicted 

Kt[  5]. 

2  A  propulsor  performance  prediction  computer  program  was  used  to  compute  the  values  of  non-dimensional 
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Kt 

K  Q 

%  Error 
Kt 

%  Error 

K  Q 

Experimental2 
MTFLOW /PBD-14.3 

0.2276 

0.2392 

0.0453 

0.0436 

5.1% 

-3.7% 

Table  5.3:  Comparison  of  propeller  4577  performance  on  Huang  Body  1  with  J  =  1.25. 

The  level  of  accuracy  for  Huang  Body  1  is  consistent  with  the  use  of  PBD/MTFLOW  as  a 
rapid,  early-stage,  design  space  exploration  tool.  Capturing  the  effective  wake  problem  in  this 
manner  allows  the  rapid  analysis  or  design  of  multiple  propellers  in  minutes  versus  hours,  as  is 
expected  with  RANS  methods. 


circulation  G  and  hydrodynamic  pitch  angle  tan  j3t  for  the  estimated  values  of  ue .  The  final  values  of  tan  were 
scaled  up  or  down  by  the  ratio  of  the  measured  value  of  Kt  to  the  computed  value  of  Kt*  The  values  of  G  and 
modified  tan/3{  were  then  used  to  compute  the  propulsor-induced  velocities  by  a  field-point  velocity  computer 
program  with  a  lifting-surface  option.  The  new  values  of  ua  were  then  used  to  compute  a  second  estimate  of 
the  effective  velocity  ue.  Three  iterations  were  sufficient  for  convergence[4]. 
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Chapter  6 


Conclusions 


6.1  Summary  of  Results 

The  coupled  PBD/MTFLOW  propeller  blade  design  and  analysis  method  is  a  viable  alterna¬ 
tive  to  axisymmetric  RANS  codes  for  solving  circumferential  mean  throughflow  problems  at  a 
significantly  reduced  computational,  as  well  as  user-preparation,  time.  The  limited  number  of 
required  starting  files  and  the  automatic  gridding  routine  of  MTSET  permits  rapid  problem 
set-up.  The  coupling  works  for  open,  ducted,  and  internal  flow  cases;  for  both  single  and 
multiblade  rows. 

6.2  Recommendations  for  Future  Work 

L  Addition  of  blade  blockage  effects  in  MTFLOW.  A  feature  of  the  tflow.xxx  file  allows 
blade  thickness  to  be  passed  into  MTFLOW  to  model  blockage,  and  thus  this  is  an  ideal 
feature  to  add  to  the  PBD2MT  program.  The  blockage  effect  becomes  important  for 
internal  multi- blade  propulsors. 

2.  Addition  of  viscous  blade  drag  in  MTFLOW.  In  the  current  coupling  method,  the  vis¬ 
cous  blade  effects  are  accounted  for  by  PBD  in  the  torque  and  thurst  calculations,  however 
the  effect  of  this  loss  on  the  fluid  velocity  is  not  passed  back  into  MTFLOW.  The  effect 
of  the  loss  is  small  on  single-blade  row  open  propellers,  but  becomes  significant  for  multi¬ 
blade  row  internal  flow  propulsors.  Accounting  for  this  loss  would  improve  the  accuarcy 
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of  the  solution  in  general,  and  would  permit  a  realistic  design  capability  for  downstream 
blade  rows  in  particular. 

3.  Fully  automated  inlet  boundary  layer  profile  modeling.  The  current  inlet  boundary  layer 
profile  method  requires  some  user  iteration  to  establish  the  entropy  loss  necessary  to 
achieve  the  desired  nominal  profile.  Altering  the  current  scaling  technique  to  a  fully 
automatic  version,  transparent  to  the  user,  would  further  reduce  the  required  user  set-up 
time  and  level  of  difficulty.  This  is  particularly  beneficial  in  that  the  modeling  of  the 
stern  section  of  a  vehicle  is  presumably  one  of  the  most  useful  aspects  of  this  coupling 
method,  and  is  thus  a  feature  which  will  presumably  be  frequently  utilized. 

4.  Improved  gridding  for  stand  alone  stern  section.  The  gridding  routine  requires  some 
manual  iteration  as  discussed  in  Section  3.4.3.  Eliminating  this  necessity  would  reduce 
the  required  set  up  time  for  stern  section  modeling. 

5.  Coupled  Propeller  Lifting  Line  (PLL) /MTFLOW.  Coupling  PLL  and  MTFLOW  would 
result  in  a  powerful  design  tool.  The  optimum  circulation  distribution  could  be  designed 
by  PLL  in  the  effective  wake  as  calculated  by  MTFLOW.  Then,  the  blade  shape  design 
could  be  accomplished  with  PBD/MTFLOW. 
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Appendix  A 


Multi-Blade  Row  Script  File 
Examples 


This  appendix  contains  example  script  files  for  a  multi-blade  row  case.  A  script  file  for  a  single 
blade  row  application  is  a  simplified  version  of  the  BatchFile  below  and  is  thus  not  included 
here. 

A.l  BatchFile 

TESTDIR=$HOME/RUNS/WJ 

ROTOR=$TESTDIR/ROTOR 

STATOR=$TESTDIR/STATOR 

bl2body=$H0ME/C0DES/MTCouple/bl2body 

velcon9=$H0ME/C0DES/VELC0N/VELC0N9/velcon9 

pbdl43=/ nf  s/ propnut /CODES/ CODES/PBD 14/PBD14 . 3/SRC/pbdl4 

pbd2mt=$H0ME/C0DES/MTCouple/pbd2mt 

mtf  lo=$HOME/CODES/Mtf  low/bin/mtf  lo 

mtsol=$HOME/CODES/Mtflow/bin/mtsol 

buildtf low=$HOME/CODES/MTCouple/buildtf low 

cd  $TESTDIR 

rm  Batch.log 
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date  >  $TESTDIR/Batch.log 
pwd  >>  $TESTDIR/Batch. log 
MAX=10 
C0UNT=1 

echo  'TITLE  =  "WJ21  PBD/MTFLOW  STATOR  COUPLING'"  »  $TESTDIR/CTQOUT . TOT 

echo  'VARIABLES  =  ' 'N'  ' ,  "Ct",  "10Cq"'  »  $TESTDIR/ CTQOUT . TOT 

echo  'ZONE  T="Ct  /  lOCq" ,  1=  '  $MAX  »  $TESTDIR/CTQOUT . TOT 

echo  'TITLE  =  "WJ21  PBD/MTFLOW  ROTOR  COUPLING'"  »  $TESTDIR/KTQOUT . TOT 

echo  'VARIABLES  =  "N"  ,  "Kt",  "10Kq"'  »  $TESTDIR/KTQOUT . TOT 

echo  'ZONE  T="Kt  /  lOKq"  ,  1=  '  $MAX  »  $TESTDIR/KTQOUT . TOT 

cd  $TESTDIR 

################################################################ 
while  test  $C0UNT  -le  $MAX 
do 

echo  '  ##  PBD  -  MTFLOW  ITERATION  '  $C0UNT  »  Batch.log 
echo  '  #######  PBD  -  MTFLOW  ITERATION  '  $C0UNT'  #######' 
################################################################ 

#  MTFLOW  TO  PBD  # 

################################################################ 

#  Create  VELJOIN.tec  from  Mtflow  results 
$bl2body  >>  Batch.log 

echo  '  ##  Done  with  b!2body  ##' 

#  Run  velcon  and  pass  req'd  files  to  ROTOR  and  STATOR  directories 
$velcon9  <  velcon. rot  >>  Batch.log 

mv  rotor.vel  $R0T0R/. 
cp  GridLower.dat  $R0T0R/. 
cp  GridUpper.dat  $R0T0R/. 
cp  GridTotal.dat  $R0T0R/. 
cp  VELJOIN.tec  $R0T0R/. 

$velcon9  <  velcon. stat  >>  Batch.log 
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mv  stator. vel  $STAT0R/. 
cp  GridLower.dat  $STAT0R/. 
cp  GridUpper.dat  $STAT0R/. 
cp  GridTotal.dat  $STAT0R/. 
cp  VELJOIN. tec  $STAT0R/ . 

################################################################ 

#  PBD  /  PBD  to  MTFLOW  # 

################################################################ 

#  Run  PBD14.3  on  Rotor  and  create  tf low. rotor 
cd  $R0T0R 

$pbdl43  <  pbd.in  >>  ../Batch. log 
echo  ’  ##  PBD  complete  on  rotor  ##* 

$pbd2mt  >>  ../Batch. log 

cp  tf low. rotor  $TESTDIR/tf low.Bl 

#cat  PBD OUT . CMV  »  PBDTOT.CMV 

#cat  PBDOUT.CMF  »  PBDTOT.CMF 

#cat  PBDOUT . SGR  »  PBDTOT.SGR 

cat  PBDOUT. KTQ  >>  ktrot.tot 

tail  -1  PBDOUT. KTQ  >  $TESTDIR/ktq. in 

read  wordl  word2  word3  word4  <  $TESTDIR/ktq. in 

echo  $C0UNT  $word3  $word4  »  $TESTDIR/KTQOUT . TOT 

#  Run  PBD14.3  on  Stator  and  create  tf low. stator 
cd  $TESTDIR/ STATOR 

$pbdl43  <  pbd.in  >>  ../Batch. log 
echo  y  ##  PBD  complete  on  stator  ##* 

$pbd2mt  >>  ../Batch. log 
cp  tf low. stator  $TESTDIR/tflow.B2 
#cat  PBDOUT. CMV  »  PBDTOT.CMV 
#cat  PBDOUT.CMF  »  PBDTOT.CMF 
#cat  PBDOUT. SGR  »  PBDTOT.SGR 
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cat  PBDGUT.KTQ  »  ctstat.tot 

tail  -1  PBDOUT . KTQ  >  $TESTDIR/ktq . in 

read  wordl  word2  word3  word4  <  $TESTDIR/ktq . in 

echo  $ COUNT  $word3  $word4  »  $TESTDIR/CTQOUT . TOT 

################################################################ 

#  MTFLOW  # 

################################################################ 

#  Combine  tflow  files  and  run  Mtflow 
cd  $TESTDIR 

$buildtflow  >>  Batch.log 
echo  *  ##  wrote  tflow. wj 21  ##* 

$mtflo  wj21exps  <  runMTFLO  >>  Batch.log 
$mtsol  wj21exps  <  runMTSOL  >>  Batch.log 
################################################################ 
C0UNT= c  expr  $C0UNT  +  lf 

#  End  of  while  loop 
done 
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A.  2  RunMTFLO 
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